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Documentation  and  Assessment  of  the  Transport  Property  Model  for  Mixtures 
Implemented  in  NIST  REFPROP  (Version  8.0) 


Justin  Chichester  and  Marcia  L.  Huber 
National  Institute  of  Standards  and  Technology* 
Boulder  CO  USA  80305-3328 


In  this  report,  we  describe  an  extended  corresponding  states  model  for 
viscosity  and  thermal  conductivity  of  mixtures  implemented  in  Version  8 of  the 
NIST  computer  program  REFPROP.  The  model  is  a modification  of  a one-fluid, 
extended  corresponding  states  (ECS)  model  for  thermal  conductivity  and  viscosity 
originally  developed  by  Ely  and  Hanley  (Ind.  Eng.  Chem.  Fundam.,  1981, 20:323- 
332).  We  apply  the  model  to  selected  mixtures  representative  of  the  fluids  available 
in  the  NIST  23  database  REFPROP  (v8),  such  as  refrigerants  and  natural  gas 
constituent  fluids,  and  present  comparisons  with  experimental  data.  Comparisons  are 
given  for  gas,  liquid,  and  supercritical  conditions. 


Keywords:  extended  corresponding  states;  natural  gas;  refrigerants;  thermal  conductivity; 
viscosity. 


* Physical  and  Chemical  Properties  Division,  Chemical  Science  and  Technology  Laboratory. 
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1.  Introduction 


The  NIST  Reference  Fluid  Thermodynamic  and  Transport  Properties  Database, 
NIST23  (REFPROP)  [1],  provides  thermophysical  properties  for  84  pure  fluids  and  for 
mixtures  of  up  to  20  components.  For  pure  fluids,  it  contains  equation  of  state  formulations 
for  the  thermodynamic  properties,  and  models  or  correlations  for  the  viscosity  and  thermal 
conductivity.  The  references  for  the  individual  pure  fluid  formulations  are  provided  in  the 
program  itself  through  the  “fluid  information”  screen.  For  mixtures,  there  is  presently  only 
one  model  available  for  the  thermal  conductivity  and  viscosity.  It  is  based  upon  the  extended 
corresponding  states  concept.  In  1981,  Ely  and  Hanley  published  a one-fluid,  extended 
corresponding  states  model  for  viscosity  [2],  and  in  1983  applied  the  model  to  thermal 
conductivity  [3].  Since  that  time,  there  have  been  several  modifications  to  the  Ely-Hanley 
method,  each  improving  upon  the  original  model  [4-10],  Variations  of  this  model  have  been 
successfully  implemented  in  two  older  NIST  databases,  NIST4  (Supertrapp)  [11]  and 
NIST14  (DDMIX)  [12],  and  most  recently  in  NIST23  (REFPROP)  [13],  In  this  report  we 
will  describe  the  formulation  used  in  the  most  recent  version  of  the  REFPROP  computer 
program  [1]  and  give  comparisons  with  experimental  data  for  the  viscosity  and  thermal 
conductivity  of  selected  mixtures  of  interest,  including  refrigerant  mixtures  and  natural  gas 
mixtures. 
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2.  Viscosity  Model 
2.1  Pure  Fluids 

In  the  extended  corresponding  states  model  used  here,  we  represent  the  viscosity  of  a 
pure  fluid  as  a sum  of  a dilute  gas  and  a residual  contribution,  and  apply  a corresponding 
states  principle  to  the  residual  contribution  only  [2], 

n(T,  p)  = rj*  (T)  + Arj(T , p)  = r/*  (T)  + Ar/0  (T0 , p0  )Fn  ( T , p), 


where  the  superscript  * denotes  a dilute  gas  value,  and  the  subscript  0 denotes  a reference 
fluid  value.  The  viscosity  of  the  reference  fluid  is  evaluated  at  a conformal  temperature  and 
density  To  and  po  given  by 


T0=T/f 


(2) 


and 

Po  = Ph- 


(3) 


The  quantities/ and  h are  called  equivalent  substance  reducing  ratios,  and  relate  the 
reference  fluid  to  the  fluid  of  interest  using  a ratio  of  critical  parameters  (denoted  by  the 
subscript  c)  and  functions  of  temperature  and  density  known  as  shape  functions  6 and  0, 

f = - ‘-0  (4) 

7,:„ 

and 


(5) 
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The  shape  factors  can  be  considered  functions  of  both  temperature  and  density.  For 
small,  nonpolar,  almost-spherical  molecules  they  are  nearly  unity  and  can  be  thought  of  as 
functions  to  compensate  for  deviations  from  a spherical  shape.  Different  approaches  have 
been  taken  to  determine  shape  factors  [4,  5,  8]  usually  depending  on  the  amount  and  quality 
of  p,V,T  information  available  for  the  fluid;  in  this  work,  we  generally  have  available 
accurate  formulations  for  the  thermodynamic  properties  of  the  pure  fluids,  either  in  terms  of 
a Helmholtz  free  energy  equation  or  a pVT  equation  of  state  (EOS),  and  we  use  a form  of  the 
“exact”  shape  factor  method  [4],  It  is  a requirement  in  this  method  to  first  determine  the 
thermodynamic  shape  factors. 


The  dilute-gas  viscosity  in  Eq.  (1)  is  given  by  the  Chapman-Enskog  theory  [14] 


n\T) 


5y]tnkBT 
16;rcrn(2’2)  ’ 


(6) 


where  the  dilute-gas  viscosity  is  r/*,  m is  the  molecular  mass,  kg  is  the  Boltzmann  constant, 
and  T is  the  absolute  temperature.  We  will  further  assume  that  a Lennard-Jones  12-6 
potential  applies,  and  use  the  Lennard-Jones  collision  diameter  for  <r.  Neufeld  et  al.  [15] 
gave  the  following  empirical  correlation  for  the  calculation  of  the  collision  integral  /22,2) 

Q<2’2)  = 1.161 45(7^*)_0 14874  + 0.524  87e”°'77320r*  + 2.161 78^243787r*  (?) 

with  the  dimensionless  temperature  T*  = k^T/s,  and  ^the  minimum  of  the  Lennard-Jones 
pair-potential  energy.  The  range  of  validity  of  this  empirical  correlation  is  0.3  < T*  < 100. 

The  factor  Fn  in  Eq.  (1)  is  found  by  using  the  expression 
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F =fU2h~2n 

where  Mis  the  molar  mass  of  the  fluid  and  M0  is  the  molar  mass  of  the  reference  fluid.  The 
model  as  developed  to  this  point  is  predictive,  and  does  not  use  any  information  on  the 
viscosity  of  the  fluid  (except  for  the  dilute-gas  piece  that  requires  Lennard-Jones  ^and  <r). 
The  functions/and  h are  found  from  thermodynamic  data  and  are  described  in  Klein  et  al. 

[9] .  In  order  to  improve  the  representation  of  the  viscosity,  an  empirical  correction  factor 
may  be  used  if  there  are  experimental  viscosity  data  available.  We  then  evaluate  Eq.  ( 1 ) at 
po,v  instead  of  po,  where  [9] 

Aw  (T,  p ) = p0(T,  p)y/(pT ),  (9) 

and  ^/is  a polynomial  in  reduced  density  px  = p! pc  of  the  form 

n 

k=Q 

(10) 

where  the  coefficients  c*  are  constants  found  from  fitting  the  experimental  viscosity  data.  As 
indicated  in  Eq.  (1),  in  order  to  evaluate  the  viscosity  of  a particular  fluid,  the  value  of  the 
residual  viscosity  of  a reference  fluid  is  required.  It  is  not  necessary  to  use  the  same 
reference  fluid  for  all  fluids.  However,  when  the  model  is  used  in  a predictive  mode,  it  is 
best  to  select  the  reference  fluid  that  is  most  similar  in  chemical  nature  to  the  fluid  of 
interest.  The  reference  fluid  should  also  have  a very  accurate  equation  of  state  and  viscosity 
surface.  When  using  pure  fluid  experimental  viscosity  to  essentially  “correct’ ’ the  viscosity, 


M 
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the  choice  of  reference  fluid  is  not  as  important,  since  an  empirical  correction  factor 
determined  from  data  is  applied,  as  in  Eqs.  (9)  through  (10). 

2.2  Mixtures 

The  extension  of  the  model  to  mixtures  is  similar  to  that  presented  for  pure  fluids,  but 
involves  an  extra  step.  First,  one  represents  the  properties  of  the  mixture  in  terms  of  a 
hypothetical  pure  fluid  (denoted  here  by  the  subscript  x)  that  is  obtained  through  the  use  of 
mixing  rules, 

/v  = /?;'Z  ZWA 

j=x  i=x  (ii) 

’ (12) 

7=1  >=\ 


and  combining  rules, 


4=VW1-'W 


(13) 


), 


(14) 


where  we  introduce  binary  interaction  parameters  and  k ^ that  can  be  determined  by 

fitting  mixture  experimental  data  when  available.  Otherwise,  they  are  set  to  zero. 
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The  iterative  method  used  to  determine  f and  hi  by  use  of  equations  of  state  is 
described  in  detail  in  the  manuscript  of  Klein  et  al.  [9],  In  that  work,  the  authors  note  that 
although  the  procedure  for  finding  the  shape  factors  is  straightforward,  it  is  complicated  in 
practice  and  the  iterative  procedure  can  fail  at  both  high  densities  and  very  low  densities.  In 
addition,  we  noted  some  failures  with  fuel  mixtures  in  the  supercritical  region  that  occurred 
with  the  use  of  the  Helmholtz  equations  of  state,  but  did  not  occur  with  a much  simpler 
equation  of  state  such  as  the  Peng-Robinson.  We  suspect  that  this  might  be  due  to  in  part  to 
the  iterations  entering  the  two-phase  region  where  the  Helmholtz  equations  can  have  very 
complicated  behavior.  To  partially  alleviate  this  problem,  we  assume  that  at  some  high 
temperature,  here  selected  as  1000  K,  the  individual  pure-fluid  shape  factors  will  have  a 
limiting  value  of 


fj 


T 

A c’7 


T. 


(15) 


h . = h. 


P, 


P 


(16) 


where  the  subscript  x denotes  the  mixture  and  j denotes  the  pure  fluid.  In  the  fluid  region 
where  the  compressibility  factor  (z=p/pRT)  of  the  fluid  mixture  is  above  0.3,  and  the 
temperature  is  above  the  critical  temperature  of  the  mixture,  we  linearly  interpolate  the 
value  of  the  f and  hj  between  its  limiting  value  in  Eqs.  (15)  through  (16)  and  the  value  at  the 
critical  temperature. 

We  assume  that  the  viscosity  of  the  mixture  obeys  a corresponding  states  principle, 
7(r,p,x)  = r/*{T,x)  + A^T^x)  = rf  (T  ,x)  + A/70(ro, Po)Fr,(T ,. p),  rlT> 
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where  the  factor  Fn  for  mixtures  is  found  with 


(18) 


n n 


1/2 


/=1  7=1 


(19) 


2 


(20) 


For  the  dilute  gas  viscosity,  binary  interaction  parameters  ky_a  and  may  be  used,  where 


(21) 


(22) 


When  there  are  sufficient  dilute-gas  viscosity  data  available,  these  parameters  are 
obtained  by  fitting  experimental  data.  Otherwise,  they  are  set  to  zero. 

3.  Thermal  Conductivity  Model 
3.1  Pure  Fluids 

We  start  with  the  procedure  of  Ely  and  Hanley  [3]  and  represent  the  thermal 
conductivity  of  a fluid  as  the  sum  of  translational  (from  collisions  between  molecules)  and 
internal  (due  to  internal  motions  of  the  molecule)  modes  of  energy  transfer, 


X{T,p)  = r\T)  + X™\T,p). 


(23) 
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The  translational  contribution  may  be  further  divided  into  a translational  dilute-gas 
contribution  (denoted  here  by  a superscript  *)  that  is  a function  only  of  temperature,  a 
residual  contribution,  and  a critical  enhancement, 

/Tans  (T , p)  = X (T)  + X (T,  p)  + Xnt  (T,  p),  (24) 

leading  to  the  following  expression  for  the  thermal  conductivity 

A(T,  p)  = Xat  (T)  + X (T)  + X (T,  p)  + X*  (T,  p).  (25) 


We  use  an  Eucken  correlation  for  the  internal  contribution 


Xnt(T)=  ^ 


M 


C — R 
p 2 


(26) 


where  Cp*  is  the  ideal-gas  heat  capacity  in  J/(mol  K),  R is  the  molar  gas  constant  [16]  (8.314 
472  J/(mol-K)),  //*  is  the  dilute-gas  viscosity  (pPa  s)  as  given  in  Eq.  (6),/int  is  set  to  1 .32- 1 O' 
3,  and  X is  in  W/(m-K).  If  sufficient  dilute-gas  thermal  conductivity  data  are  available,  fmX  is 
fit  to  a polynomial  in  temperature, 

/int  = a0  + aiT-  (22) 

For  the  dilute-gas  translational  contribution  (in  W/(m-K))  we  use 

X(T)=  15x10  R'1  , (28) 

4 M 

where  the  dilute  gas  viscosity  ;/*  is  from  Eq.  (6).  The  residual  contribution  is  found  using 
extended  corresponding  states: 

X(T,p)  = A’0(T„,pl))Fx,  (29) 

with 


9 


FA=fU2h~ 2/3 


Mo 

M 


(30) 


In  order  to  improve  the  representation  of  the  thermal  conductivity,  an  empirical 
correction  factor  may  be  used  if  there  are  experimental  thermal  conductivity  data  available. 
We  then  evaluate  Eq.  (29)  at  po^  instead  of  po , where  [8] 

PoAT’P)  = Po(T’P)x(Pr\  (31) 


and  / is  a polynomial  in  reduced  density  pT  = plpQ  of  the  form 

z(p,)  = t.b>P^  (32) 

k=  o 


where  the  coefficients  /y  are  found  from  fitting  the  experimental  thermal  conductivity  data. 

The  critical  contribution  is  computed  using  a simplified  crossover  model  developed 
by  Olchowy  and  Sengers  [17], 

pC  R0kBT  , N 

a (T,p)=  - f - -(n-n„),  (33) 

OTTTjg 


where  the  heat  capacity  at  constant  pressure,  Cp(T,p),  is  obtained  from  the  equation  of  state, 
Rq  = 1 .03  is  a universal  constant  [18],  and  the  viscosity, ;/( T,p),  is  obtained  from  the  method 
described  earlier.  The  crossover  functions  Q and  Qo  are  determined  by 


n = 


n 


( C -C  ^ 

p ' 


c 


p y 


arctan  (qd%)  + ^r(qd£) 


(34) 
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(35) 


n„  = - 

n 


1 - exp 


-1 


(^r1  + 


(q&Pc 

p 


The  heat  capacity  at  constant  volume,  Cv(T,p),  is  obtained  from  the  equation  of  state,  and  the 
correlation  length  £ is  given  by 


PcP 

v/ y 

dp{T,p) 

Tr  dp(TR,p) 

lrA2J 

dp 

T T dp 

T _ 

(36) 


The  partial  derivative  of  density  with  respect  to  pressure  dp/dp\r,  is  evaluated  with  the 
equation  of  state  at  the  system  temperature  T and  a reference  temperature,  7r.  The  reference 
temperature  is  a value  where  the  critical  enhancement  is  assumed  to  be  negligible:  7r  = 
1.5TC.  The  exponents  y = 1.239  and  v = 0.63  are  universal  constants  [18],  The  critical 
amplitudes  r and  Co  are  system-dependent  and  can  be  determined  by  the  asymptotic  behavior 
of  the  equation  of  state  in  the  critical  region.  For  most  fluids,  we  set  these  to  the  fixed  values 
r=  0.0496  and  = 1.94  x 10'10  m.  The  thermal  conductivity  at  the  critical  point  itself  is 
infinite.  The  cutoff  wavenumber  qa  (or  alternatively,  its  inverse  q^)  is  system-dependent 
and  determined  from  regression  of  experimental  data  when  available  or  estimated  [19]. 

3.2  Mixtures 

The  extension  of  the  model  to  mixtures  is  similar  to  that  presented  for  viscosity, 
A(T,  p,  x ) = X (T,  x)  + Amt  (T,x)  + AX  (T,  p)  + Xnt  (T,  p,  x)  (37) 

and  only  the  residual  contribution  is  treated  with  corresponding  states. 

A(T,  p,  x)  = X (T,  x)  + A,nt  (T,  x)  + AA0  (T0 , p0  )FA  (T,p)  + /T  ( T,  p,  x) . (38) 
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The  internal  and  translational  dilute-gas  contributions  for  the  mixture  are  found  with  the 


empirical  mixing  rule, 


« x;(/t'nt(r)  + /i ](T)) 

'C,  (T,*)  + 2mlx  <T,  x)  = X 1 

i=i 


(39) 


with 


(1  - )(l  + (A^y  / A/,)’'4)2 

(m  + Mj/M,)}12 


(40) 


All  quantities  are  evaluated  at  the  mixture  temperature  T,  and  the  dilute-gas  viscosity  is 
found  with  Eq.  (6).  The  parameter  klhx  is  an  empirical  binary  interaction  parameter  for  the 
dilute  gas  region  determined  by  fitting  experimental  data. 


The  residual  contribution  for  the  mixture  requires  calculation  of  a mixture  Fx  with 


Fx=fl‘2K™g"\ 


Ml/2  " - 

Jx  .T  /•= 1 j= 1 


-1/2 


gl/2  = 

/-1/2;.  4/3 


1/2 


1 Igi+Vg, 


(kF, 


J 


m 

(42) 


with  gi  found  with 


g,  = MC 


W’Pj) 


-4/3 


J J 


(43) 


Similar  to  the  procedure  for  the  calculation  of  the  viscosity,  we  use  the  combining  rules 
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(44) 


(45) 


where  we  introduce  binary  interaction  parameters  and  kijthx  that  can  be  determined  by 
fitting  experimental  data  when  available.  Otherwise,  they  are  set  to  zero. 


As  observed  for  pure  fluids,  there  also  is  an  enhancement  of  the  thermal  conductivity 
observed  near  the  critical  point;  however,  theory  for  mixtures  is  not  as  well  developed  as  for 
pure  fluids.  In  REFPROP  8 we  make  a simple  approximation  where  we  obtain  the 
parameters  necessary  for  the  simplified  Olchowy-Sengers  model  (Eqs.  33  through  36,  above) 
by  taking  a mole  fraction  average, 


Tk  =1.52>,rCi  (46) 

/=! 


</d=2>d,  <47) 

i= 1 


r = I*,r, 

m (48) 

=!>,£•  *49) 

i= 1 

The  other  parameters  y , v are  universal  and  are  unchanged  from  their  pure-fluid  values.  The 
thermodynamic  properties  are  evaluated  at  the  mixture  T and  p,  not  at  scaled  values,  and  the 
mixture  Tc,pc  and  pc  are  used.  In  addition,  at  the  critical  point  the  thermal  conductivity  of  a 
pure  fluid  is  infinite  [20];  for  a mixture,  theory  and  experiment  suggest  that  a finite  value  is 
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achieved  [21].  Therefore  we  have  imposed  a somewhat  arbitrary  upper  limit  (100  % of  the 
background  value)  on  the  magnitude  of  the  enhancement  possible  in  a mixture. 


4.  Results 


We  tested  the  model  on  a limited  set  of  mixture  data  containing  both  refrigerant 
mixtures  and  hydrocarbon  mixtures.  When  possible,  we  obtained  binary  interaction 
parameters  to  improve  the  representation  of  the  data.  The  parameters  were  obtained  by  least 
squares  regression  of  experimental  data.  The  resulting  values  for  the  binary  interaction 
parameters  for  both  the  dilute  and  residual  viscosity  are  given  in  Table  1 below.  These 
binary  interaction  parameters  are  contained  in  the  file  HMX.BNC  that  is  distributed  with  the 
REFPROP  computer  program.*  In  addition,  we  provide  a list  of  the  Lennard- Jones  12-6 
parameters  for  pure  fluids  in  REFPROP8  in  Table  2.  References  for  these  values  are 
provided,  however,  often  it  was  necessary  to  estimate  the  values.  The  estimation  methods  of 
Chung  et  al.  [22]  and  a corresponding-states  estimation  method  [7]  were  used. 

Table  3 gives  the  deviations  of  the  model  compared  with  experimental  data  for 
selected  binary  systems.  We  use  the  following  definitions  for  average  absolute  deviation 
(AAD),  bias,  and  root-mean-square  (RMS)  deviation: 


(50) 


(51) 


and 
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(52) 


RMS 2 = 


100 


n 

I 

/=l 


.2\ 


7, 


exp 


-1 


100  ^ 


( 'N'N 


z 


n, 


" ,•= 


exp 


-1 


J) 


In  general,  the  model  represents  the  viscosity  of  the  mixtures  to  within  5 % to  10  %.  One 
notable  exception  is  the  data  for  the  C02/decane  system  measured  by  Cullick  and  Mathis 
[23].  The  model  does  not  represent  this  system  well  even  with  the  use  of  interaction 
parameters.  This  is  probably  due  to  the  large  size  differences  between  the  molecules.  It  has 
been  noted  previously  [24]  that  one-fluid  corresponding  states  types  of  models  have 
difficulties  representing  systems  where  the  molecules  differ  greatly  in  size.  The  trend  of 
increasing  deviations  as  the  size  difference  between  the  molecules  increases  is  also 
demonstrated  by  the  data  of  Aucejo  et  al.  [25]  in  the  systems  of  pentane  with  increasingly 
larger  alkanes,  with  the  largest  deviations  observed  for  the  pentane/dodecane  system.  Further 
work  in  this  area  is  needed. 

There  also  are  large  discrepancies  between  the  model  predictions  and  the  viscosity 
data  of  Heide  [26]  for  the  binary  mixture  of  R 143a  and  R125.  This  system  has  also  been 
studied  by  Laesecke  [27]  and  is  represented  well  by  the  model,  so  we  suspect  that  there  may 
have  been  an  experimental  problem  with  the  Heide  [26]  data  for  this  system. 
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Table  la.  Values  of  interaction  parameters  for  viscosity  for  selected  refrigerant, 
cryogen  and  hydrogen  mixtures. 


System 

kij,a 

kij.e 

kij.fii 

kij.h'l 

Source  of 
data 

R134a/R32 

0.0437 

-0.4678 

-0.0645 

-0.0305 

[281  [291  [301 

R 125/propane 

-0.21 

0.3776 

0.0 

0.0 

[30] 

R22/propane 

-0.1284 

0.3331 

0.0 

0.0 

[311 

R1 15/propane 

-0.2744 

0.414 

0.0 

0.0 

[311 

R3  2/propane 

0.2857 

-0.4346 

0.0 

0.0 

[28] 

R12/R22 

-0.0213 

0.0698 

0.0 

0.0 

[31] 

R13/R22 

-0.0332 

0.1389 

0.0 

0.0 

[3H 

R13B1/R22 

-0.2505 

0.7028 

0.0 

0.0 

[311 

R14/R22 

0.0523 

-0.1194 

0.0 

0.0 

[311 

R 14/CO 

-0.0421 

-0.0638 

0.0 

0.0 

[32] 

R22/R152a 

-0.1382 

0.4355 

0.0 

0.0 

[311 

R32/R124 

0.0 

0.0 

-0.0813 

-0.1158 

[291 

R125/R134a 

-0.0045 

0.0038 

-0.0829 

-0.0106 

[281  [291  [331 

hydrogen/methane 

-0.2744 

-0.1738 

0.0 

0.0 

[22,  341 

hydrogen/nitrogen 

-0.6225 

0.1388 

0.0 

0.0 

[34] 

hydrogen/C02 

-0.9218 

0.8531 

0.0 

0.0 

[34J 

hydrogen/oxygen 

-0.7661 

0.3683 

0.0 

0.0 

[341 

hydrogen/CO 

-0.7599 

0.8536 

0.0 

0.0 

[341 

hydrogen/helium 

-0.0480 

-1.9733 

0.0 

0.0 

[34] 

hydrogen/argon 

-1.0538 

0.8177 

0.0 

0.0 

[34] 

hydrogen/neon 

-0.6105 

-0.4449 

0.0 

0.0 

[341 

hydrogen/ammonia 

-0.5275 

0.8419 

0.0 

0.0 

[34] 

hydrogen/ethylene 

-0.7131 

0.8778 

0.0 

0.0 

[34] 

hydrogen/propylene 

-0.5744 

0.7662 

0.0 

0.0 

[34] 

hydrogen/cis-butene 

-0.7844 

0.8518 

0.0 

0.0 

[341 

helium/argon 

-0.187 

-2 

0.0 

0.0 

T351 

krypton/CO 

-0.08889 

-0.05653 

0.0 

0.0 

[32] 

16 


Table  lb.  Values  of  interaction  parameters  for  viscosity  for  selected  natural  gas 
component  and  hydrocarbon  mixtures. 


System 

kij.cr 

kij,£ 

kijfn 

kij.hti 

Source  of 
data 

methane/CO 

-0.0165 

-0.2493 

0.0 

0.0 

[32] 

methane/C02 

-0.00677 

-0.3342 

0.0 

0.0 

T361 

methane/ethane 

0.0331 

-0.2549 

0.0 

0.0 

[371 

methane/propane 

-0.0182 

-0.0263 

0.0 

0.0 

[371 

methane/butane 

0.0030 

-0.1292 

0.0 

0.0 

[361 

methane/nitrogen 

-0.3594 

0.9076 

-0.0612 

0.0398 

[38] 

ethane/propane 

0.0557 

-0.2792 

0.0 

0.0 

[371 

ethane/butane 

0.0350 

-0.1858 

0.0 

0.0 

[37] 

ethane/hydrogen 

-0.4544 

0.3877 

0.0 

0.0 

r 34i 

propane/butane 

0.0825 

-0.3471 

0.0 

0.0 

[371 

propane/hydrogen 

-0.5926 

0.7304 

0.0 

0.0 

[34] 

pentane/heptane 

0.0 

0.0 

-0.0892 

-0.1474 

[251 

pentane/dodecane 

0.0 

0.0 

-0.1188 

-0.4911 

[251 

pentane/decane 

0.0 

0.0 

-0.0965 

-0.2775 

[25] 

pentane/nonane 

0.0 

0.0 

-0.1195 

-0.1610 

T251 

hexane/octane 

0.0 

0.0 

-0.0519 

-0.0010 

[25] 

hexane/dodecane 

0.0 

0.0 

-0.0168 

-0.2999 

[251 

hexane/decane 

0.0 

0.0 

-0.0550 

-0.1310 

[251 

hexane/nonane 

0.0 

0.0 

-0.0263 

-0.0738 

[25] 

hexane/cyclohexane 

0.0 

0.0 

-0.12 

0.2243 

[391 

heptane/dodecane 

0.0 

0.0 

-0.0716 

-0.1600 

[25] 

heptane/decane 

0.0 

0.0 

-0.0374 

-0.0677 

[251 

heptane/nonane 

0.0 

0.0 

-0.0704 

-0.0117 

[251 

octane/dodecane 

0.0 

0.0 

-0.0802 

-0.1160 

[251 

octane/decane 

0.0 

0.0 

-0.0299 

-0.0488 

[251 

nonane/dodecane 

0.0 

0.0 

-0.0634 

-0.0803 

[251 

nonane/decane 

0.0 

0.0 

-0.0671 

-0.0241 

[251 

decane/dodecane 

0.0 

0.0 

-0.0499 

-0.0614 

[251 
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Table  2.  Lennard-Jones  12-6  parameters  for  selected  pure  fluids  in  REFPROP  v8 


Fluid 

£^b(K) 

o'(nm) 

Ref. 

Fluid 

s^b(K) 

rr(nm) 

Ref. 

Rll 

363.6 

0.5447 

Estimated 

ammonia 

386 

0.2957 

[40] 

R12 

297.24 

0.5186 

Estimated 

argon 

143.2 

0.335 

[41] 

R13 

204.00 

0.4971 

[42] 

butane 

280.51 

0.573 

[43] 

R14 

164.44 

0.4543 

[44] 

CO 

91.7 

0.369 

[12] 

R32 

289.65 

0.4098 

Fit  to  [45] 

C02 

251.196 

0.3751 

[46] 

R115 

201.90 

0.5876 

[47] 

cyclohexane 

297.1 

0.6182 

[43] 

R116 

226.16 

0.5249 

Estimated 

H20.D20 

809.1 

0.2641 

[43] 

R124 

275.80 

0.5501 

[48] 

decane 

490.51 

0.686 

[49] 

R125 

249.00 

0.5190 

[50] 

dimethyl  ether 

395 

0.4307 

[43] 

R 1 4 1 b 

370.44 

0.5493 

Estimated 

ethane 

245 

0.4362 

[51] 

R142b 

278.20 

0.5362 

[52] 

ethanol 

362.6 

0.453 

[43] 

R218 

266.35 

0.5800 

Estimated 

ethylene 

224.7 

0.4163 

[43] 

R227ea 

289.34 

0.5746 

Estimated 

H2S 

301.1 

0.36237 

[43] 

R236ea 

318.33 

0.5604 

Estimated 

helium 

10.22 

0.2551 

[43] 

R236fa 

307.24 

0.5644 

Estimated 

heptane 

400 

0.64947 

[12] 

R245ca 

345.44 

0.5505 

Estimated 

hexane 

399.3 

0.5949 

[43] 

R245fa 

329.72 

0.5529 

Estimated 

hydrogen 

59.7 

0.2827 

[43] 

RC318 

299.76 

0.5947 

Estimated 

isobutene 

332 

0.5026 

Estimated 

R113 

376.035 

0.6019 

Estimated 

isohexane 

395.2 

0.5799 

Estimated 

Krypton 

178.9 

0.3655 

[43] 

isopentane 

341.06 

0.56232 

[12] 

Methane 

174 

0.36652 

[53] 

isobutane 

307.55 

0.46445 

[54] 

Methanol 

481.8 

0.3626 

[43] 

nitrogen 

98.94 

0.3656 

[41] 

N20 

232.4 

0.3828 

[43] 

nonane 

472.127 

0.66383 

[49] 

Neon 

32.8 

0.282 

[43] 

oxygen 

118.5 

0.3428 

[41] 

neopentane 

191 

0.644 

[55] 

propane 

263.88 

0.49748 

[56] 

Octane 

452.09 

0.63617 

[49] 

R23 

243.91 

0.4278 

[57] 

Pentane 

341.1 

0.5784 

[12] 

R41 

244.88 

0.4123 

Estimated 

Propylene 

298.9 

0.4678 

[43] 

R114 

323.26 

0.5770 

Estimated 

R22 

284.724 

0.4666 

fit  to  [58] 
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Table  3.  Summary  of  results  for  the  viscosity  of  binary  mixtures  for  selected  systems. 
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Table  4.  Summary  of  results  for  the  viscosity  of  multicomponent  mixtures. 


Max 

3.88 

-2.26 

-1.62 

10.71 

-4.31 

-2.47 

-AM 

-6.69 

-8.50 

-6.75 

RMS 

1.12 

0.54 

0.30 

r- 

rn 

04 
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1.14 
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960 

1.07 

1.09 

BIAS 

1.28 
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-1.20 
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CO 
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-5.05 

-6.18 

-4.30 

AAD 

rn 

0.71 
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1.42 

2.38 
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1 91 
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6.18 

4.30 

npts 

157 

224 

224 
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34 
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Laesecke[28] 

Schley[65] 

Schley[65] 

VC 

fc 
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Carr[66] 
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Assael[67] 

Langelandsvik[68] 
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3 
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na 

968  0 

^r 

oo 

o 

0.735 

0.731 

0.956 

0.848 

0.902 

O 

O 

oo 

o 

0.922 

p range,  kPa 

221-2853 

95-20073 

95-20227 

101-69154 

101-66051 

101-58364 

242-14040 

4460-25110 

13310-25310 

4950-25240 

T range,  K 

246-340 

260-320 

260-320 

298-339 

299-394 

303-398 

241-455 

263-303 

263-304 

263-304 

System 

refrigerant  mixture 
(R32/R125/R134A) 

natural  gas  H 

natural  gas  L 

natural  gas  1 

natural  gas  2 

natural  gas  3 

natural  gas 

natural  gas  1 

natural  gas  2 

natural  gas  3 
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We  also  made  comparisons  of  the  viscosity  prediction  of  the  REFPROP  model  for 
multicomponent  mixtures.  Comparisons  with  a ternary  refrigerant  mixture  R32/R125/R134A 
measured  by  Laesecke  et  al.  [28]  are  presented  in  Table  4 and  Figures  la-c.  In  the  Figures, 
the  percent  deviation  is  defined  as  1 00*(calculated  value-experimental)/experimental.  Two 
different  compositions  0.33/0.33/0.34  and  0.3/0. 1/0.6  mole  fraction  of  R32/R125/R134a 
were  measured  from  246  to  340  K at  pressures  up  to  2853  kPa  with  average  absolute 
deviations  of  less  than  two  percent.  The  estimated  uncertainty  of  this  data  is  given  by  the 
author  as  2.4  %.  We  also  made  comparisons  with  several  natural  gas  mixtures.  Figures  2a-c 
show  the  deviations  between  the  REFPROP  model  and  experimental  data  for  several 
different  natural  gas  mixtures  as  a function  of  density,  pressure  and  temperature.  Schley  et 
al.  [65]  measured  two  different  natural  gases  at  pressures  to  20  MPa.  The  "L"  gas  is  a low 
calorific  gas  with  a relatively  high  content  (9.75  mole  %)  of  nitrogen;  the  high  calorific  "FT 
natural  gas  contains  only  1.5  % nitrogen.  Both  samples  are  represented  well  by  the 
REFPROP  model.  We  note  that  a binary  interaction  parameter  for  methane/nitrogen 
viscosity  was  added  to  the  REFPROP  program  after  the  initial  release  of  Version  8 in  April 
2007;  this  parameter  significantly  improved  the  results  for  the  high-nitrogen  natural  gas. 
Carr  [66]  measured  three  natural  gas  samples.  Sample  1 is  a high-ethane  natural  gas  (73.5 
mole  % methane,  25.7  % ethane,  0.6  % N2  and  0.2  % propane),  Sample  2 is  a high-nitrogen 
natural  gas  (73. 1 % methane,  15.8  % N2,  6. 1 % ethane,  3.4%  propane,  0.8  % helium,  0.6  % 
n-butane  and  0.2  % isobutane)  obtained  from  a transmission  line,  and  Sample  3 is  a 
simulated  low-ethane  natural  gas  (95.6  % methane,  3.6  % ethane,  0.5  % propane  and  0.3  % 
N2).  With  the  exception  of  the  Sample  1 gas  at  pressures  above  20  MPa,  the  model 
represents  the  data  well.  The  Assael  et  al.  [67]  natural  gas  is  a 5 component  mixture  of 
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methane,  ethane,  propane,  CO2  and  nitrogen  ( 84.8,  8.4,  0.5,  5.6,  and  0.66  mol  % 
respectively)  at  pressures  up  to  14  MPa  and  is  represented  with  an  average  absolute 
deviation  of  less  than  two  percent,  but  shows  systematic  deviations  with  the  largest 
deviations  at  the  lowest  temperatures.  Systematic  deviations  are  even  more  pronounced  in 
the  data  from  Langelandsvik  et  al.  [68]  who  measured  the  viscosity  of  natural  gas  samples 
from  three  different  locations  in  the  North  Sea.  The  primary  constituents  in  gas  1 are: 
methane  90.2  %,  ethane  6.3  %,  propane  0.8  % and  CO2  1.8  %;  in  gas  2:  methane  80.0  %, 
ethane  9.3  %,  propane,  5.0  %,  CO2  2.2  %;  and  in  gas  3:  methane  92.2  %,  ethane  4.4  %, 
propane  0.5  %,  CO2  1 %.  REFPROP  shows  systematic  negative  deviations  for  all  the 
Langelandsvik  et  al.  [68]  samples  with  the  deviations  increasing  as  the  pressure  increases  at 
the  lowest  temperatures.  It  was  noted  in  the  Langelandsvik  et  al.  [68]  manuscript  that  there 
are  unexplained  differences  between  the  measurements  of  Schley  et  al.  [65]  and  the  work  of 
Langelandsvik  et  al.  [68].  The  REFPROP  program  is  in  better  agreement  with  the 
measurements  of  Schley  et  al.  [65], 


The  thermal  conductivity  model  was  tested  on  a limited  set  of  thermal  conductivity 
mixture  data  containing  both  refrigerant  mixtures  and  hydrocarbon  mixtures.  When  possible, 
we  obtained  binary  interaction  parameters  to  improve  the  representation  of  the  data.  The 
resulting  values  for  the  binary  interaction  parameters  for  both  the  dilute  and  residual  thermal 
conductivity  are  given  in  Table  5 below. 
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Table  5.  Values  of  interaction  parameters  for  thermal  conductivity  for  selected 
systems. 


System 

kn,a 

knM 

Data  source 

methane/C02 

-0.15 

0.0 

0.0 

169] [70] 

methane/CO 

-0.07 

0.0 

0.0 

[71] 

methane/propane 

-0.2 

0.0 

0.0 

[721 

methane/ethane 

-0.12 

0.0 

0.0 

[731 

R22/R142b 

-0.76 

0.0 

0.0 

[74]  [75] 

R32/propane 

0.18 

-0.08 

0.09 

[76] 

R32/R134a 

0.08 

-0.09 

0.1 

[76] 

R32/R125 

0.0 

-0.15 

0.16 

[77,  78]  [79] 

R134a/propane 

0.07 

-0.04 

0.20 

[761 

Comparisons  for  both  binary  and  multicomponent  mixtures  were  made.  The  results 
are  summarized  in  Tables  6-8  below.  Most  systems  that  do  not  contain  data  near  the  critical 
region  have  deviations  of  less  than  5 % with  the  notable  exception  of  the  binary  system 
R32/R125  measured  by  Kim  et  al.  [80],  which  shows  deviations  over  10  %.  Since  this  same 
binary  system  was  measured  by  other  researchers  [77-79]  and  the  model  agrees  well  with  the 
experimental  data,  there  must  be  an  unexplained  problem  with  the  data  of  Kim  et  al.  [80] 
for  this  binary  system.  Perkins  et  al.  [76]  made  extensive  measurements  on  binary  mixtures 
containing  R32,  R125,  R134a,  and  propane.  The  measurements  covered  the  gas,  liquid  and 
some  near-critical  points.  Figures  3a-d  show  the  deviations  for  Perkins  et  al.  [76] 
measurements  as  a function  of  density.  The  REFPROP  model  represents  the  data  well  except 
for  points  where  the  critical  enhancement  is  significant;  in  these  cases  the  model 
significantly  underestimates  the  thermal  conductivity. 

Roder  and  Friend  [73]  made  extensive  measurements  on  the  thermal  conductivity  of  a 
binary  mixture  of  methane  and  ethane  for  three  different  concentrations  (0.7,  0.5,  and  0.35 
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mol  fraction  methane)  over  a range  of  conditions  that  included  regions  where  the 
enhancement  is  significant.  Sakonidou  et  al.  [2 1 ] also  measured  the  methane/ethane  system 
in  the  critical  region.  Figures  4a-d  show  the  percent  deviations  of  the  data  from  the  model 
with  and  without  the  critical  enhancement.  The  figure  demonstrates  the  improvement  seen 
when  the  enhancement  term  is  included;  however,  it  also  demonstrates  that  the  present 
model  is  not  adequate,  and  further  work  in  this  area  is  needed. 

5.  Conclusions 

We  describe  a modified  corresponding  states  model  based  on  an  extended 
corresponding  states  model  originally  developed  by  Ely  and  Hanley  [2,  3]  for  the  viscosity 
and  thermal  conductivity  of  mixtures.  The  model  has  been  implemented  in  the  computer 
program  REFPROP  [1],  version  8,  available  from  NIST.  Comparisons  with  selected  mixture 
data  for  viscosity  and  thermal  conductivity  are  presented.  The  model  represents  the  viscosity 
and  thermal  conductivity  of  mixtures  to  within  5 to  1 0 %,  except  for  mixtures  containing 
large  size  differences,  where  the  deviations  are  larger.  A model  for  the  thermal  conductivity 
of  mixtures  in  the  critical  region  is  also  implemented,  and  although  it  provides  some 
enhancement  to  the  thermal  conductivity,  comparisons  with  experimental  data  indicate  that 
the  model  significantly  underpredicts  the  thermal  conductivity  in  the  critical  region,  and 
further  research  in  this  area  is  recommended. 


We  thank  our  NIST  colleagues  Dr.  Richard  Perkins  and  Dr.  Daniel  Friend  for 
their  helpful  comments.  JC  acknowledges  a Professional  Research  Experiences  Program 
(PREP)  undergraduate  fellowship  at  NIST. 
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Table  6.  Summary  of  results  for  the  thermal  conductivity  of  binary  mixtures  for  selected  systems. 
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Table  7.  Summary  of  results  for  the  thermal  conductivity  of  multicomponent  mixtures  for  selected  systems. 
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Figure  la.  Percent  deviation  of  viscosity  as  a function  of  temperature  for  a mixture  of 
refrigerants  R32,  R125  and  R134a. 
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Figure  lb.  Percent  deviation  of  viscosity  as  a function  of  density  for  a mixture  of 
refrigerants  R32,  R125  and  R134a. 
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Figure  lc.  Percent  deviation  of  viscosity  as  a function  of  pressure  for  a mixture  of 
refrigerants  R32,  R125  and  R134a. 
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Figure  2a.  Percent  deviation  of  viscosity  as  a function  of  density  for  selected  natural  gas 
mixtures. 
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Figure  2b.  Percent  deviation  of  viscosity  as  a function  of  pressure  for  selected  natural  gas 
mixtures. 


33 


% Deviation 


15 


10 


5 


0 


-5 


-10  L 
200 


♦ Assael  et  ai.  2001 

■ Carr  1953  Natural  Gas  1 
A Carr  1953  Natural  Gas  2 
X Carr  1953  Natural  Gas  3 
x Schley  et  al.  2004  Natural  Gas  H 

• Schley  et  al.  2004  Natural  Gas  L 

+ Langelandsvik  et  al.  2007  Natural  Gas  1 
- Langelandsvik  et  al.  2007  Natural  Gas  2 
Langelandsvik  et  al.  2007  Natural  Gas  3 


A 


250  300  350  400  450  500 

Temperature,  K 


Figure  2c.  Percent  deviation  of  viscosity  as  a function  of  temperature  for  selected  natural 
gas  mixtures. 
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Figure  3a.  Percent  deviation  of  thermal  conductivity  as  a function  of  density  for  binary 
mixtures  of  propane  and  R32. 
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Figure  3b.  Percent  deviation  of  thermal  conductivity  as  a function  of  density  for  binary 
mixtures  of  R 134a  and  R32. 
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Figure  3c.  Percent  deviation  of  thermal  conductivity  as  a function  of  density  for  binary 
mixtures  of  propane  and  R134a. 
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Figure  3d.  Percent  deviation  of  thermal  conductivity  as  a function  of  density  for  binary 
mixtures  of  R 125  and  R134a. 
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Figure  4a.  Percent  deviation  of  thermal  conductivity  as  a function  of  density  for  binary 
mixtures  of  methane  and  ethane  near  the  critical  region  with  a critical  enhancement  term 
in  the  model. 
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Figure  4b.  Percent  deviation  of  thermal  conductivity  as  a function  of  density  for  binary 
mixtures  of  methane  and  ethane  near  the  critical  region  without  a critical  enhancement 
term  in  the  model. 
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Figure  4c.  Percent  deviation  of  thermal  conductivity  as  a function  of  temperature  for 
binary  mixtures  of  methane  and  ethane  near  the  critical  region  with  a critical 
enhancement  term  in  the  model. 
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Figure  4d.  Percent  deviation  of  thermal  conductivity  as  a function  of  temperature  for 
binary  mixtures  of  methane  and  ethane  near  the  critical  region  without  a critical 
enhancement  term  in  the  model. 
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